#include <vtkPoints.h>
#include <vtkPointData.h>
#include <vtkPolyData.h>
#include <vtkSmartPointer.h>
#include <vtkFloatArray.h>

int main(int, char *[])
{
  vtkSmartPointer<vtkPoints> points =
    vtkSmartPointer<vtkPoints>::New();
  points->InsertNextPoint(1,1,1);
  points->InsertNextPoint(2,2,2);
  points->InsertNextPoint(3,3,3);

  vtkSmartPointer<vtkPolyData> polydata =
    vtkSmartPointer<vtkPolyData>::New();
  polydata->SetPoints(points);

  vtkSmartPointer<vtkFloatArray> floatArray =
    vtkSmartPointer<vtkFloatArray>::New();
  floatArray->SetNumberOfValues(3);
  floatArray->SetNumberOfComponents(1);
  floatArray->SetName("FloatArray");
  for(vtkIdType i = 0; i < 3; i++)
    {
    floatArray->SetValue(i,2);
    }
  polydata->GetPointData()->AddArray(floatArray);

  vtkSmartPointer<vtkIntArray> intArray =
    vtkSmartPointer<vtkIntArray>::New();
  intArray->SetNumberOfValues(3);
  intArray->SetNumberOfComponents(1);
  intArray->SetName("IntArray");
  for(vtkIdType i = 0; i < 3; i++)
    {
    intArray->SetValue(i,2);
    }

  polydata->GetPointData()->AddArray(intArray);

  for(vtkIdType i = 0; i < 3; i++)
    {
    double p[3];
    polydata->GetPoint(i,p);
    vtkFloatArray* pointsFloatArray = vtkFloatArray::SafeDownCast(polydata->GetPointData()->GetArray("FloatArray"));
    vtkIntArray* pointsIntArray = vtkIntArray::SafeDownCast(polydata->GetPointData()->GetArray("IntArray"));
    std::cout << "Point " << i << " : " << p[0] << " " << p[1] << " " << p[2] << " "
              << pointsFloatArray->GetValue(i) << " " << pointsIntArray->GetValue(i) << std::endl;
    }

  polydata->GetPointData()->NullPoint(1);
  polydata->Modified();

  for(vtkIdType i = 0; i < 3; i++)
    {
    double p[3];
    polydata->GetPoint(i,p);
    vtkFloatArray* pointsFloatArray = vtkFloatArray::SafeDownCast(polydata->GetPointData()->GetArray("FloatArray"));
    vtkIntArray* pointsIntArray = vtkIntArray::SafeDownCast(polydata->GetPointData()->GetArray("IntArray"));
    std::cout << "Point " << i << " : " << p[0] << " " << p[1] << " " << p[2] << " "
              << pointsFloatArray->GetValue(i) << " " << pointsIntArray->GetValue(i) << std::endl;
  
    }

  return EXIT_SUCCESS;
}
